A template bank to search for gravitational waves from inspiralling compact binaries: 

II. Phenomenological model 
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Matched filtering is used to search for gravitational waves emitted by inspiralling compact binaries 
in data from ground-based interferometers. One of the key aspects of the detection process is the 
deployment of a set of templates, also called a template bank, to cover the astrophysically interesting 
region of the parameter space. In a companion paper, we described the template-bank algorithm 
used in the analysis of LIGO data to search for signals from non-spinning binaries made of neutron 
star and /or stellar-mass black holes; this template bank is based upon physical template families. 

In this paper, we describe the phenomenological template bank that was used to search for 
gravitational waves from non-spinning black hole binaries (from stellar mass formation) in the 
second, third and fourth LIGO science runs. We briefly explain the design of the bank, whose 
templates are based on a phenomenological detection template family. We show that this template 
bank gives matches greater than 95% with the physical template families that are expected to be 
captured by the phenomenological templates. 

PACS numbers: 02.70.-c, 07.05.Kf, 95.85.Sz, 97.80.-d 



I. MOTIVATION 

The Laser Interferometer Gravitational- Wave Obser- 
vatory (LIGO) detectors [l[ have reached their design 
sensitivity curves. The fifth science run (S5) began in 
November 2005 and should be completed by the end of 
2007, with the goal of acquiring a year's worth of data 
in coincidence between the three LIGO interferometers. 
Each successive LIGO science run has witnessed improve- 
ment from both experimental and data analyst's point of 
view. On the experimental side, better stationarity of 
the data and detector sensitivities closer to design sen- 
sitivity curve were achieved. On the data analysis side, 
the search pipeline was tuned, and new techniques were 
developed to reduce the background rate while keeping 
detection efficiencies high. 

Among the sources of gravitational waves that ground- 
based detectors are sensitive to, inspiralling compact bi- 
naries are among the most promising. Several searches 
for inspiralling compact binaries in the LIGO data have 
been pursued: primordial black holes (PBH) 0, bi- 
nary neutron stars (BNS) [H, 0, JBL, an d intermediate 
mass binary black holes (BBH) [g, Q. These searches 
used matched filtering technique, which is the most effec- 
tive and commonly used method to search for inspiralling 
compact binaries. 

Matched filtering computes cross correlation between 
the detector output and a template waveform. If the tem- 
plate waveform is identical to the signal then the method 
is optimal, in the sense of signal-to-noise ratio (SNR). 
However, in general, the template waveforms differ from 
the signals. Indeed, modelizations can only approximate 
the exact solution of the two-body problem. In addi- 
tion, template waveforms are constructed with a subset of 
the signal's parameters (e.g., the two component masses 
whereas eccentricity and spins effects may be neglected) . 
In this work, we consider only the case of non-spinning 



waveforms so that the signals are entirely defined by 4 
parameters: 2 mass components, the time of arrival and 
the initial phase. Because the signal's parameters are un- 
known, the detector output must be cross correlated with 
a set of template waveforms, which is called a template 
bank. While the spacing between templates can be de- 
creased most certainly, and this is the insurance of a SNR 
close to optimality, it also increases the size of the tem- 
plate bank (i.e., the computational cost). The distance 
between templates is governed by the trade-off between 
computational cost and loss in detection rate; therefore, 
template bank placement is a key aspect of the detection 
process. 

In a companion paper 0] , we proposed a template bank 
with a minimum match of 95%. We assume that both 
template and signal were based on the same physical tem- 
plate family (precisely, the stationary phase approxima- 
tion with the phase described at 2PN order |g, Hj] ) . We 
have shown that the template bank could be used, effec- 
tively, to search for BNS, PBH, black hole - neutron star 
(BHNS), and BBH systems. In this paper, we consider 
the case of BBH systems only. 

There is a wide variety of techniques used to describe 
the gravitational flux and energy generated during the 
late stage of the inspiralling phase (e.g., see [HI). How- 
ever, they lead to various physical template families, and 
overlaps between them are not necessarily high. In the 
case of heaviest systems, post-Newtonian (PN) expan- 
sion [l(| begins to fail as the characteristic velocity v/c 
is not close to zero anymore (e.g., see [ll|). In addi- 
tion, even though heavier BBH systems are accessible 
each time the detector sensitivity improves in the low 
frequency range, BBH waveforms remain short in the 
LIGO band. For instance, during the second science run 
(S2) @, the lower cut-off frequency was set to 100 Hz, 
which restricted the total mass of the search to be be- 
low 40 Mq, and the longest expected signal to last about 
0.60 s. 
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There exist several template families, and there is no 
reason to select one in particular. A solution may be to 
filter the detector output with a set of template banks, 
each of them associated with a different physical template 
families. We have shown in [l7j that a unique template 
bank placement could be used effectively with several 
template families. However, we investigate only 4 dif- 
ferent families at 2PN order. A different template bank 
might be necessary for other template families. More 
importantly, the number of template families could be 
large, and computational cost not manageable. Instead 
of searching for BBH signals using several physical tem- 
plate families, a single detection template family (DTF) 
was proposed by Buonanno, Chen, and Vallisneri [11| 
(BCV) with the goal of embedding the different phys- 
ical approximations all into a single phenomenological 
model. This detection template family is known as BCV 
template family and has been used to search for non- 
spinning BBH signals in LIGO data [3, @| ■ 

In this paper, we do not intend to compare a search 
that uses BCV templates and a search based upon phys- 
ical template family. Our main goal is to describe the 
BCV template bank that was developed and used to 
search for stellar-mass BBH signals in the second_[S2), 
third (S3), and fourth (S4) LIGO science runs @,[|. In 
section HH we briefly discuss the template parameters and 
the filtering process related to BCV templates. In sec- 
tion IIIIl we describe the BCV template bank, and the 
the spacing between templates. In section IIV| we test 
and validate the proposed template bank with exhaustive 
simulated injections. Finally, in scctionfVl we summarize 
the results. 



II. THE BCV TEMPLATE FAMILY 

The detection template family that was proposed 
in [ill ] is built directly from the Fourier transform [l2j 
of gravitational-wave signals by writing the amplitude 
and phase as polynomials in the gravitational-wave fre- 
quency law that appear in the stationary phase approx- 
imation [T3|. In the frequency domain, the BCV tem- 
plates arc defined to be 



h(f)=A(f)e^, 



(2.1) 



where 



A(f) = r V6 (l - «/ 2/3 ) (/cut - /) , (2.2) 

N 

i>(f) = 2irft + 4>o + r 5/3 Y;^f k/3 ■ ( 2 - 3 ) 

The parameters to and <po are the standard time of arrival 
and initial phase of the gravitational wave signal. The 
parameter a is a shape parameter introduced to capture 
post-Newtonian amplitude corrections. Because various 
models predict different terminating frequencies, an end- 
ing cut-off frequency / cu t is introduced. In the ampli- 



tude expression, the waveform is multiplied by a Hcav- 
iside step function, 8 (/ cut — /). In the right hand side 
of equation 12.31 we use only two parameters ipo and ip3, 
which suffices to obtain a high match with most of the 
PN models [ll|. The symbol ^3 here is the same as the 
symbol ip3/2 i n [HI- 

The iftk parameters arc the phase parameters of the 
phenomenological waveform, which cannot be directly 
linked to the physical mass parameters; the BCV tem- 
plates are made for detection, not for parameter esti- 
mation. Nevertheless, a good approximation (for low 
masses) of the chirp mass M. is given by 
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In section IIVD1 we investigate the range of validity of 
this relation. 

The filtering of a data set using BCV templates is not 
as trivial as the one that uses physical template fami- 
lies. Indeed, the BCV filtering implies a search in six 
dimensions (-00, 4*3, ct, 4>o, to, and /cut)- The SNR can 
be analytically maximized over a, <j)o and to, which re- 
duces the search to three dimensions only, the maxi- 
mization over a and (f>o is not. In order to perform 
the filtering and the maximization over a and <po, we 

need to construct orthonormal basis vectors \ hk > 

I J fc=l,..,4 

for the 4-dimensional linear subspace of templates with 
4>o £ [0,27r) and a £ (—00,00), and we want the basis 

vectors to satisfy (hi hj^j = Sij (see the Appendix for 

details). The SNR before maximization is given by 

p = Xi cos ui cos 4>o + X2 sin to cos <f>o (2-5) 
+ X3 cos a; sin (/>o + x\ sin a; sin 0o , 

where = /s /kjY and s is the data to be filtered 1 . The 

parameter w is a function of a (sec cquation lA8l aiid lA9 1. 

The SNR p can be maximized over cf> and oj(a). In [6], 
the maximization is done over the two new parameters 
A = u> + </>o and B = lu — </>q. The maximized SNR 
(independent of a and (fio), is denoted pu, and is given 
by 



Pu 




(2.6) 



where Vk are function of Xi (see Appendix and equations 
lAl3 | lAlil and lATi]) . 

The SNR provided in cquation l2.6l is the unconstrained 
SNR that is independent of any constraint on the range 
of the parameter cP F = a/ 2 / 3 (again, here the index U 



The expression of the SNR shows that the expected rate of false 
alarm follows a chi-square distribution with 4 degrees of freedom 
instead of 2 in the case of physical template families 
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represents the unconstrained case, and we shall use C for 
the constrained case). Yet, in [Tl(, the authors suggested 
that the parameter ap should be restricted to the range 
[0, 1]. Indeed when a v F > 1, the amplitude in eauation l2.2l 
becomes negative, which corresponds to unphysical wave- 
forms. Moreover, when ap < 0, the amplitude factor can 
substantially deviate from the predictions of PN theory. 

In S2 many accidental triggers were found with 
ckp > 1, and the calculation of the SNR was uncon- 



strained (as in cquation l2.6[) leading to a high false alarm 
rate, which was decreased, a posteriori, by removing all 
triggers for which ctp > 1 (without decreasing the de- 
tection efficiency). Nevertheless, triggers that verified 
a l p < were kept because the false alarm rate did not 
decrease significantly when this selection was applied. 

In S3 and S4, the search for BBH systems deployed 
a filtering that takes value into account, by using a 



maximization of equation 12 . 51 that leads to a constrained 
SNR denoted pc- The expression for the constrained 
SNR depends now on the value of ctp. We have ap = ap 
if < oF F < 1, and no constraint is applied (i.e., pc = 
Pij). However, if oF F < or dtp > 1, then a constrained 
SNR is used so that the final atp parameter is and 
1, respectively, and pc < Pu- The expressions of the 
constrained SNR are provided in the Appendix. 

Let us make the point clear. In the contrained-SNR 
case, as explained above, ap* = ap if < < 1 only. 
Otherwise ap takes only two values (0 or 1) but we know 
ap (since it is the condition to apply the constraint or 
not). 

It is worth noticing that for the study that follows, 
we always use a constrained SNR but using an uncon- 
strained SNR should not significantly change the results 
of our simulations and/or template bank placement. In- 
deed, simulated injections are generated with physical 
template families for which we do not expect olf to be un- 
physical (i.e., outside [0—1]), as we shall see in sectionlVl 
Using a constrained SNR has an important impact when 
dealing with real analysis, where most of the acciden- 
tal triggers have ap between (—00,00) (and therefore 
in (— 00, 0[ U ]l,oo) as well, where pc < Pu)- Conse- 
quently, in general, for a given threshold A, the SNR of 
accidental triggers have pc < pu, and the rate is there- 
fore lower with respect to a search with unconstrained 
SNR. The number of triggers that needs to be stored is 
lower by an order of magnitude. Nevertheless, the fi- 
nal rate of triggers between the two methods may be 
equivalent because of a posteriori cuts on aj? when an 
unconstrained SNR is used as in lal. 



III. BCV TEMPLATE BANK DESIGN 



Template bank placement has been investigated in sev- 
eral papers @, EH, EH, El, E3, EH m the context of physical 
template families. We refer the reader to the established 
literature in this subject area. 



A. Metric Computation in Plane 



In the case of BCV templates, the mismatch metric 
9ij [HI is known (see Appendix), and is constant over 
the entire ipo - ^ parameter space. Nonetheless, the met- 
ric components are strongly related to the lower cut-off 
frequency of the search, which affects the moments used 
to calculate the metric (see equation IB6j) . The moment 
computation also depends on the a parameter, as dis- 
cussed later. For now, let us suppose that the moments 
are fixed. 

Because the metric is constant, the placement of tem- 
plates on the ifjQ-Tpz parameter space is straightforward. 
In the first search for BBH signals Q , the template place- 
ment used a square lattice, and templates were placed 
parallel to the ipo axes. In S3 and S4 BBH searches, an 
optimal placement was used (hexagonal lattice), which 
reduced the requested computing resources (and trigger 
rate) by 30% with respect to S2. In this paper, we only 
consider tests related to the hexagonal lattice case. In S3 
and S4 BBH searches, we placed the templates parallel to 
the first eigenvector rather than parallel to the ?/>o axis. 

The target waveforms are BBH systems for which the 
lowest component mass is set to 3 Mq and the high- 
est component mass is defined by the detector lower 
cut-off frequency (up to 80 Mq in S4). Simulations 
show that to detect such target waveforms, the range 
of phenomenological parameters should be set to V'o € 
[10000, 550000] Hz 5/3 and V > 3 e [-5000, -10] Hz 2/3 . 
As explained in section HlID| if we search for BBH sys- 
tems only, a significant fraction of those templates are 
not needed and can be removed from the template bank. 



B. as-dependence 



The moments used to estimate the metric components 
strongly depend on the parameter a. We refer to this 
parameter as as to differentiate from the a parameter 
(or equivalcntly from ap) that is used in the filtering 
process. As shown in figure [TJ the number of templates 
changes significantly when cxb varies. There is a drop in 
the number of templates around 0.3 = 10~ 2 . We want 
to minimize the template bank size but we also need to 
consider the efficiency of the bank as defined in 0, E3 , 
and choose cxb appropriately. Indeed, we expect the ef- 
ficiency of the template bank to be also affected by this 
parameter. We performed simulated injections so as to 
test the efficiency of the template bank for various values 
of ag. Results are summarized in figure [5] for three typi- 
cal values of cxb- Because efficiencies are very similar, we 
decided to use anag parameter such that the number of 
templates is close to a minimum, that is 10~ 2 . In all the 
following simulations and LIGO searches, as = 10 -2 . 
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FIG. 1: The size of the template bank function of the o.b pa- 
rameter. The two curves show the template bank size versus 
as parameter for two values of sampling frequency. The two 
curves show the same pattern with a drop around olb = 1CP 2 , 
where template bank sizes are twice as low as compared to 
Ob = or olb = 2.5 x 10~ 2 . This evolution of the template 
bank size is directly linked to the moment computation (see 
equation IB6[) . where the parameter olb is used. Efficiency of 
a template bank is not strongly related to this parameter (see 
figure[2]) so we choose a value that corresponds to the smallest 
bank size. In real analysis, we use 2048 Hz and for simplicity 
a value of olb = 10 -2 was chosen. In this example, we used 
the same simulation parameters as in section llV Al 



C. Template Bank using ending frequency layers 

Starting from each template that is placed in the ipo-Tps 
plane, we need to lay templates along the third dimen- 
sion, which is the ending cut-off frequency / cut of the tem- 
plate. Because the mismatch is first order in A/ cut fill ], 
it cannot be described by a metric. 

Using an exact formula, [Tl| proposed to lay tem- 
plates with different / cu t values between / m i n and / en d = 
/Nyquist that depend on the region searched for. We pop- 
ulate the / cu t dimension as follows. First, we estimate 
the frequency of the last stable orbit which we refer to as 
/min, and the frequency at the light ring which we refer to 
as / max . Between / min and / max , we place N cut layers of 
templates with the ending frequency chosen at equal dis- 
tance between / m ; n and / ma x- The frequency at the last 
stable orbit and light ring are defined in terms of the total 
mass (/ min = 1/(A/tt6 3 / 2 ), / max = 1(Mtt3 3 / 2 )). The to- 
tal mass is computed for each template using an empirical 
expression similar to equation 12.41 M « —ip3/(S2TT 2 ipo). 
This expression is an approximation. It underestimates 
the total mass for low mass range, however, it is suitable 
for the wide range of mass that we are interested in: the 
final template bank gives high match with the various 
physical template families, as shown in section ITVl In all 
our simulations and searches, we set the minimal match 
(MM) [l5| to 95%, so there is no guarantee that the rela- 
tion between M and ^0,3 is suitable for minimal matches 
far from 95%. 
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FIG. 2: Template bank efficiencies versus aj parameter. The 
as parameter does not significantly affect the matches. Most 
of them are above the minimal match of 95%, and more im- 
portantly the three distributions are close to each other. In 
this example, we used the same simulation parameters as in 
section HV Al and EOB injections. 
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FIG. 3: Template bank efficiencies versus number of layers, 
A cut , in the / cut dimension. With the current template bank 
design, N cu t does not affect the matches significantly. The 
cumulative distribution of matches over 10,000 simulations 
shows only small differences between 3 and 20 layers. Even 
the results obtained with 1 layer are not that far from A cut = 
3. In our real analysis and simulations, we used Acut = 3. In 
this example, we used the same simulation parameters as in 
section HV Al and EOB simulated injections. 



D. Polygon Fit 

The boundaries of the template bank are defined by 
the ranges of the parameters ipi and the span of the 
cut-off frequency / cut in such a way that BBH systems 
with component mass as low as 3 M Q arc detectable. 
The ipi ranges provided in section IIII Al cover a squared 
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area that is actually too wide: a significant fraction of 
the templates are not targeting the BBH systems we are 
searching for. Therefore, in order to reduce the template 
bank size and optimize our searches, we introduce an ex- 
tra procedure that selects the pertinent templates only. 
This procedure is known as a polygon fit and works as 
follows. First, we create a BCV template bank with the 
range of tpo and ip3 parameters as large as possible, and 
for our purpose, as quoted in section UlI Al This choice 
of ranges allows us to not only detect BBH systems but 
also BHNS systems. Since, we focus on the BBH systems 
only, we perform many BBH simulated injections and fil- 
ter them with the template bank that has been created. 
For each injection, we keep the ipo and -03 parameters 
of the template that gives the best match. We gather 
all the final pairs of 'ijjQ, ip3 parameters, and superpose 
them on top of the original template bank. It appears 
that only about a third of the templates are required to 
detect BBH systems with a high match. This sub-set 
of templates can be used to define a polygon area that 
enclose all of them. The resulting polygon area defines 
the boundaries of our new template bank and results in a 
template bank three times smaller than the original one. 

In figurelU we show such a template bank that is within 
the boundary of a polygon constructed with our sim- 
ulated injections. The coordinates of this polygon are 
chosen empirically. For safety, the boundaries are chosen 
loosely, therefore the template bank has also the ability 
to detect non-spinning BHNS. It is worth noticing that 
with this template bank designed to detect BBH many 
BHNS systems are found with a match greater than the 
requested MM (See section HVC|) . 




FIG. 4: Example of parameter space and template bank place- 
ment. This plot shows a projection of the templates onto the 
ipo /4>3 plane. Simulations and equation 12.41 gives an estima- 
tion of the mapping between the phenomenological parame- 
ters and the chirp mass of the simulated injections. Low mass 
systems such as a (3,3)Mq are in the RHS, high mass sys- 
tems lie on the LHS and asymmetric systems in the bottom 
left corner. In this example, we used the same parameters as 
in section llVBI 



IV. SIMULATIONS 

In the following simulations, we fix the sampling fre- 
quency to 2048 Hz, as = 10~ 2 , N cut = 3, the ipo and if>3 
ranges are provided in section IIII Al and a polygon fit as 
in figure 2] is used. The simulated injections are based on 
several physical template families that are labelled EOB, 
PadeTl, TaylorTl, and TaylorT3 0, HSJp, HI HI with 
the phase expressed at 2PN order (see [17| f° r more de- 
tails). The population of simulated injections has a uni- 
form total mass. Although this choice is not based on any 
astronomical observation, it is convenient to estimate the 
efficiency of our template banks. We use a noise model 
that mimics the design sensitivity curve of initial LIGO 
(see [HG3), and tne minimal match is MM = 95%. We 
performed 2 simulations that are closely related to the 
third and fourth LIGO science run's BBH searches @ . 

A. Example 1 

The first set of simulations uses a lower cut-off fre- 
quency of 70 Hz, as in S3 BBH search The maximal 
total mass of the simulated injections is set to 40 M© 
and therefore the largest component mass to 37 M@. 
The template bank has 531 templates. The results are 
summarized in figure [5] which shows the efficiency of the 
template bank versus the total mass. There are a few in- 
jections found with a match as low as 93% for total mass 
M < 6.5 M Q . Closer inspection shows that several issues 
are linked to this feature. First, we used a sampling fre- 
quency of 2048 Hz, which reduces the template bank size 
by « 50% as compared to a sampling of 4096 Hz. Second, 
we set olb = 10~ 2 , which reduces the template bank size 
by sa 50% as compared to «b =0. Finally, the number 
of layers, A^ cu t, is limited to 3. Therefore, this tuning sig- 
nificantly reduced the template bank size with the cost 
of losing about only 1 to 2% SNR for a small fraction 
of the parameter space considered. From M = 6.5 M© 
to about M = 20 M Q , matches are above 95%. In the 
high mass range, a large fraction of the simulated in- 
jections are found below the minimal match (but larger 
than 90%): 20% in the case of TaylorTl, TaylorT3, and 
PadeTl models, and only 0.1% in the case of EOB injec- 
tions. This effect is expected because the lower cut-off 
frequency is high, and therefore many of the high mass 
systems considered are very short (i.e., less than 100 ms). 
Because the final frequency of the EOB signals goes up to 
the light ring, the matches are larger than in the case of 
TaylorTl, TaylorT3, and PadeTl approximants, whose 
last frequencies stop at the last stable orbit. 

B. Example 2 

The second set of simulations uses a lower cut-off fre- 
quency of 50 Hz, as in S4 BBH search The maximal 
total mass of the simulated injections is set to 80 M© and 
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FIG. 5: Distribution of the efficiency versus the total mass. 
The simulation consisted of N s — 40, 000 injections. The 
lower cut-off frequency for the injections and the BCV tem- 
plates was set to 70 Hz, as in S3 BBH search. 



therefore the largest component mass to 77 Mq. The 
template bank has 1609 templates. The results are sum- 
marized in figure [Sj Up to M « 40 M Q , most of the 
injected simulations are recovered with matches above 
95%. However, a small fraction is found with matches 
below 95%, which represent 0.1% of the EOB, PadeTl, 
and TaylorTl injections, and 3% of the TaylorT3 injec- 
tions. In the high mass region (up to 60 Mq), 20% of the 
injections are below the required minimal match for the 
TaylorTl, TaylorT3, and PadeTl injections, and only 
0.5% of the EOB injections. If we consider injections 
with total mass from 60 to 80 Mq, almost 10% of EOB 
are below the minimal match (but above 92%). As for 
other models, matches drop quickly towards zero down 
to 40%, which is due to shorter and shorter duration of 
the injected waveforms. 



FIG. 6: Distribution of the efficiency versus the total mass. 
The simulation consisted in N s = 40, 000 injections. The 
lower cut-off frequency of the injections and the BCV tem- 
plates is 50 Hz, as in S4 BBH search. For convenience (same 
scale as in figure [SJ), we do not show simulated injections with 
total mass above 60 Mq and matches below 90% (see text for 
details) . 



60% of the BHNS injections are recovered with the match 
larger than 95%, 77% with the match larger than 90%, 
and 98% with the match greater than 50%. Therefore, 
using the same bank as in S3 and S4 searches, whose 
boundaries resulting from the polygon fit were deliber- 
ately chosen to be slightly wider than necessary, we can 
detect a significant fraction of the BHNS systems. It is 
also clear from the figure that the lightest systems have 
a very low match. This was expected since the tem- 
plate bank aimed at detecting systems whose total mass 
is greater than 6Mq, as defined by the maximum of the 
-0o range. 



C. Example 3 

As stated in section UlIDl although the template bank 
is designed to target BBH systems, it has the ability to 
detect some BHNS systems as well. The goal of this third 
simulation is to demonstrate that indeed many BHNS 
systems are detectable with a high match by using the 
template designed to search for BBH systems in S3 and 
S4 data sets. 

The parameters used are exactly the same as in the 
second example. The maximal total mass of the simu- 
lated injections is set to 80 Mq, the largest component 
to 79 Mq, and the lowest component mass is set to 1 Mq. 
We impose the systems to be BHNS only (the mass of the 
neutron star is less than 3 Mq , and the mass of the black 
hole is larger than 3 Mq). The template bank is identical 
to the second simulation (1609 templates). The results 
are summarized in figure where we plot matches as a 
function of the two component masses. We found that 



We performed a second test where the polygon fit is 
not applied anymore. The template bank is then much 
larger with 4635 templates but we found that 78% of the 
BHNS injections are recovered with a match larger than 
95%, 94% with a match larger than 90%, and 98% with 
a match greater than 50%. The size of such a template 
bank is comparable to a template bank that uses phys- 
ical template families (e.g., with the same parameters 
as above, a hexagonal placement for physical template 
families |17| that covers a parameter space from 1 to 80 
solar mass has about 3000 templates if we exclude the 
templates for which both component mass are below 3 

Mq). 

The events which are found with a low match (say, 60% 
or lower) correspond to low mass systems where the neu- 
tron star's mass is less than 2.5 Mq and the BH's mass 
less than 7 Mq which can be taken care of by increasing 
the range of ipo. 



7 





2.8 




2.6 








2.4 


CO 

- 


2.2 


i 






2 


v. 




a 


1.8 


c 






1.6 




1.4 




1.2 



' ■ 1 ' t 


I 

3 


1 n filial*! « ''' 

■Kr ft L " r\ 


cJ ; 

'" ' 1 \ 
it i * 1 



























10 20 30 40 50 60 70 

Black hole mass M n 



-0.9 
0.8 
■0.7 
0.6 
0.5 
0.4 
lo.3 
lo.2 
'o.1 



FIG. 7: Matches between the BCV template bank used to 
search for BBH systems and BHNS injections. The simulation 
consisted in N 3 — 100, 000 injections. We found that 77% of 
the BHNS injections have a match larger than 90%. See the 
text for more details. 
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FIG. 9: Errors in chirp mass estimation corresponding to sim- 
ulations made in section IIVB I Errors increase significantly for 
A4 greater than 8 solar mass with errors larger than 50%. 




a second peak in the range [1, 1.1], which correspond to 
about 15% of the injections; it corresponds to total mass 
above 60 Mq. 

In figure [51 we plot the errors on the chirp mass (i.e., 
(Mi — A4 e )/M.i, where i stands for injected, and e for 
estimated). We used equation 12.41 to estimate the chirp 
mass. The errors are within 10% for BBH systems when 
the chirp mass is below as 8 Mq. However, errors increase 
significantly when M. > 8 Mq because (i) parameter es- 
timation of high mass BBH systems is intrinsically weak, 
even for physical template families and (ii) BCV tem- 
plates are known to be detection template families that 
are not suitable for parameter estimation. 



FIG. 8: Distribution of the parameter corresponding to 
simulations made in section IIVBI Most of the found sim- 
ulated injections have an a% value between zero and unity. 
However, a significant fraction are distributed around aj. — 1. 
Those triggers correspond to M > 60 Mq, for which wave- 
forms cannot be differentiated from a transient noise (short 
duration) . 



D. Discussions 



In this section, we use the results of section IIVBI to 
check (i) the range of validity of equation l2.41 which gives 
an estimation of the chirp mass, and (ii) the regime of 
constrained SNR (i.e., the value of of). 

Although we use a constrained SNR, we kept track 
of the value of o^ F before the maximization. We plot 
the distribution of oF F in figure [5] About 83% of the 
injections were found with a dtp value in the range ]0,1[. 
Therefore, as stated in section [H] the results obtained 
with the constrained SNR are very similar to what we 
would have obtained if we had used the unconstrained 
SNR. The distribution has a first peak around 0.7 and 



V. CONCLUSIONS 

The BCV template bank that we described in this pa- 
per was used to search for BBH systems in the S2, S3 
and S4 LIGO data sets. We described the significant im- 
provements that were made between the S2 search and 
the S3/S4 searches: olb tuning, hexagonal lattice, and 
polygon fit. These improvements reduce the template 
bank size by an order of magnitude, while keeping the effi- 
ciency higher than 95% for most of the BBH systems con- 
sidered. Consequently, despite reducing the lower cut-off 
frequency from 100Hz to 50Hz between the second and 
the fourth science runs, the template bank size remained 
similar. 

A principal motivation for the construction of a detec- 
tion template bank was to use a single template bank 
instead of several physical template families. The tem- 
plate bank size is therefore an important aspect of a BCV 
search, and we have shown in this work how the num- 
ber of templates can be optimized to search for BBH. 
Remarkably, the same template bank has a high match 
with a wide range of BHNS. 
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More importantly, the BCV template bank was de- 
signed to search for BBH systems in the context of the 
S2 LIGO search. That is. for a lower cut-off frequency 
of 100 Hz for which most of the target waveforms are 
short duration waveforms. However, LIGO detectors im- 
proved and are still improving at low frequency, making 
the waveforms longer. The advantage of using a BCV 
template to search for systems as low as 3 M is no longer 
evident, especially considering the absence of a well de- 
fined x 2 test for phenomenological templates. Therefore 
if it were to be used, the author thinks that a BCV tem- 
plate bank should be used to search for a mass range 
starting at a higher value, such as 10 or 20 Mq. 



Acknowledgments 

We would like to acknowledge many useful discussions 
with members of the LIGO Scientific Collaboration, in 
particular of the LSC- Virgo Compact Binary Coalescence 
working group, which were critical in the formulation of 
the results described in this paper. This work has been 
supported in part by Particle Physics and Astronomy 
Research Council, UK, grant PP/B500731. This paper 
has LIGO Document Number LIGO-P070089-01-Z. 



APPENDIX A: FILTERING, a-MAXIMIZATION, 
AND CONSTRAINED SNR 



We define the inner product as follows 

hi(m(f) 



hi 



4Re 



df 



(Al) 



/o S h (f) 

where Sh (/) is the one-sided noise power spectral density. 



1. Filtering 

The BCV templates in the frequency domain are de- 
fined by equation 12.11 The amplitude part of a BCV 
template A(f) can be decomposed into linear combina- 
tions of /~ 7 / 6 and / _1 / 2 . These expressions can be used 
to construct an orthonormal basis {hk}k=i.A- We want 
the basis vectors to satisfy 



(A2) 



First, we construct two real functions Ai(f) and A2(f) 
that satisfy (A^A,) = Sij. Then, we define hi^if) = 
Ai,a(/)e <t *"* ,1 l ft3,4(/) = iAi a {fy^-^ which will 



give 



hi\h: 



Sij, and the desired basis, {h^}. ip is 



the phase of the signal, as defined in equation 12.31 and 
(j>o the initial phase that we want to maximize. We can 
choose the following basis functions 







a n 




7" 


-7/6 " 


Mf)_ 




0,21 0,22 




./" 


-1/2 



where the normalization factor are given by 



7—1/2 

on — i 7 ^ 3 , 02i 




-1/2 



-1/2 



(A3) 

(A4) 
(A5) 

-*7/3 J 

and the integrals Ik arc defined by 

/"/cut f— k 

4 = 4 /. sefj* (A6) 

The normalized template can be parametrized using 
the orbital phase </>q and an angle oj 

h(6,oj;f) = (A7) 

hi (/) cos cj cos <j>o + /i2(/)sinu;cos(/)o 

+ hs(f) cos a; sin (/)o + hi(f) sin a; sin ^>o 

where w is related to a by (see (f|) 

on& , . „ N 

tanw = , (A8) 

022 + a 2 ia 

which can be inverted to get a 

a 22 tanw 

a = . (A9) 

a2i tancj + an 

It follows that for any given signal s, the overlap is 



P 



(A10) 



Xi COS 0J COS <po + 3?2 sin 0J COS 0o 

£3 cos oj sin cf>o + X4 sin oj sin c/)q , 



where Xi = (^s hi). We can then maximized over oj (i.e., 

a), and 0o without any constraint on the a parameter, 
which leads to the unconstrained SNR given by 



Pu 




where 



V 
Vi 
V 2 



x\ 
x\ 



x 2 



1 

x 2 



X 2 , 

T 2 



2(xiX 2 + £324). 



(All) 



(A12) 
(A13) 
(A14) 



The values of oj and 4>o that maximize pu are provided 
in Q as function of the Xi . We reformulate 
the Vi and found this simple expression: 



^max usm g 



Vn 

tan2w max = — 
Vi 



(A15) 



It is then straightforward to obtain a max using equation 



91 and 



a 



max y2/3 
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constrained and unconstrained SNRs 



Starting from equation IA101 we can derive a con- 
strained SNR pc that depends upon the value of the 
parameter ap. Therefore, we need to maximize equa- 
tion [XlO] over the parameter 4>o only. This maximization 
gives 



1 



p(uj) = max — F i/Vb + V\ cos 2lo + ~Vi sin 2uj . (A16) 

<" V2 

If qj? < 0, we want to use a SNR calculation for which 
lo = 0, which means = a = 0. Therefore, the con- 
strained SNR is 



Finally, using the relation V\ cos 2lo + V2 sin 2lu = 
\/Vi + V 2 2 cos (2lj — 6*), where tanfl = V 2 /Vi, we can re- 
write equation |A16] in the general case where < otp < 1 
by imposing 2lu = 9, which gives 



PC 




(A20) 



that is an identical expresion as in equation IA11I (i.e., 
Pc = pu)- So, equation IA15I is also valid, and off = otp. 
More details on this derivation can be found in 1241. 



p° c = —^V + V 1 



(A17) 



If olf > 1, we want to use a SNR calculation for which 



which means that a. 



1 (i.e., a = fj t /3 ). 
Using equation IA8| the angle u> = cj m ax is then a maxi- 
mum given by 



arctan ■ 



,-2/3 

a lJcut 



The constrained SNR is then given by 



Pc 



V2 



cos 2cj max + V2 sin 2uj„ 



(A18) 



(A19) 



APPENDIX B: METRIC 



We can derive an expression for the match between 
two BCV templates (described by equation 12.31 12,ll and 
12. 2p . First, we consider templates with the same ampli- 
tude function (i.e., the same a and / cut parameter). The 
overlap (h(ip , -0 3 ), h(ipo + A?/'o, "03 + ^^3)} between tem- 
plates with close values of tpo and ip3 can be described (to 
second order in Ai[>o and A-03) by the mismatch metric 

9ij & 



(/i(V>o,V , 3),MV ) o + bipods + M's)) = 1 - MiMj- 

i,i=0,3 



(Bl) 



r 



The metric coefficients can be evaluated analyti- 
cally [ll|, and are given by 



where uq = 5/3 and n 3 = 2/3, and 



1 



9ij 



[Mi - MjMz l M 2 ] 



where the Mn),„(3) are the matrices defined by 



M 



(i) 



M 



(2) 



M 



(3) 



J(2n ) J(n + n 3 ) 
J(n + n 3 ) J{2n 3 ) 

J (no) Jim) 
J(no-l) J(n 3 -1) 

J(0) J(-l) 
J(-l) J(-2) 



(B2) 

(B3) 
(B4) 
(B5) 



J(n) = 



kK/)l 2 1 



4f 



Sh(f) 
that the 



Let us emphasize the fact 
(h(tpo,ip3),h(ip ■ 

variant in the ipQ-ip 3 plane, so the metric g^ 



'If 



(B6) 
mismatch 



"03 + Aip 3 )} is translationally in- 
is constant 

everywhere since J{n) is independent of V'chV's parame- 
ters. 
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